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The in-medium physics of heavy quarkonium is an ideal proving ground for our ability 
to connect knowledge about the fundamental laws of physics to phenomenological pre- 
dictions. One possible route to take is to attempt a description of heavy quark bound 
states at finite temperature through a Schrodinger equation with an instantaneous po- 
tential. Here we review recent progress in devising a comprehensive approach to define 
such a potential from first principles QCD and extract its, in general complex, values 
from non-perturbative lattice QCD simulations. Based on the theory of open quantum 
systems we will show how to interpret the role of the imaginary part in terms of spatial 
decoherence by introducing the concept of a stochastic potential. Shortcomings as well 
as possible paths for improvement are discussed. 
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1. Motivation 

The investigation of the quark gluon plasma created in relativistic heavy-ion cohi- 
sions rehes on a fruitful interplay between theory and experiment. Hence the avail- 
ability of observables that are both experimentally accessible as well as theoretically 
amenable is of importance. The first unambiguous sign for the creation of a decon- 
fined and strongly interacting state of matter, reported by the RHIC experiments 
PHENIX^ and STAR^, e.g. relied on the occurrence of collective flo^^ ^ and the 
quenching of jets 133 in Au -h Au collisions at y^^NN = 200GeV. 

Early on, Matsui and Satz'' proposed another type of hard probe, the bound 
states of heavy quarks. From analogies with the electromagnetic plasma, they con- 
cluded that a Debye screened potential above the deconfinement temperature would 
not be able to support the formation of a cc bound state. This in turn would lead 
to a suppression of the measured abundance of heavy quarkonia in the presence of 
a quark-gluon plasma. 

While heavy quarkonium suppression has indeed been confirmec£Mi3III by sev- 
eral experiments, it turns out, that the physics of charmonium in relativistic heavy- 
ion collisions is rather involved. Cold nuclear matter effects, as well as final state 
effects, such as recombination and feed-down mask the actual influence of the ther- 
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mal medium on the stability of the bound state. 

In the case of 66, which became measurable with adequate statistics, once the 
LHC commenced its runs at y^^NN = 2.76TeV, the connection is much cleaner. Pro- 
duced in the hard partonic stage at the beginning of the collision, bottomonium does 
not suffer from significant recombination effects nor is there appreciable feed-down 
from tt. The dimuon spectra measured by the CMS collaboratioiP^ furthermore pro- 
vided clear evidence for the suppression of the excited states in Pb + Pb collisions 
compared to their abundances in p + p. Under these conditions heavy quarkonium 
becomes an ideal candidate to investigate the physics of melting bound states, which 
will be our focus in the rest of this brief review. 

Ever since Matsui and Satz, the goal for theory has to be to put their phenomeno- 
logical arguments onto a solid first-principles basis. In the past many attempts relied 
on a time independent notion of the melting process, introducing purely real model 
potentials and investigating their ability to support a bound state of Q and Q. 
Popular candidates were the color singlet free energies F^(r)P^, the color singlet 
internal energies even linear combinations of both quantitieJ^^. While 

these quantities can be readily extracted from lattice QCD correlation functions at 
a single time step /3 = r, a direct connection to a Schroedinger equation derived 
from first principles QCD is unlikely to exist^^. 

In a situation where the heavy quarkonium is produced inside the quark-gluon 
plasma, e.g. through recombination, questions about the possibility of bound state 
formation can lead to insight. Bottomonium however is created before the QGP 
thermalizes and we need to understand its real-time evolution inside the medium 
to eventually grasp the physics of its suppression. 

An essential step into this direction was the discovery by Laine et.al.^that the 
heavy quark potential at first non-trivial order in resummed perturbation theory not 
only shows Debye screening but also features an imaginary part. Its appearance was 
subsequently attributed to the phenomenon of Landau dampin^I^. First corrections 
to these results in the context of effective field theory were presented i]JI^ introducing 
further contributions to both real and imaginary part. These findings remind us that 
a static notion of a well defined bound state above the deconfinement temperature 
becomes devoid of meaning and urge us to devise a fully dynamic approach to the 
stability and melting of heavy quarkonium instead. 

Let us begin by asking how a non-relativistic Schroedinger equation with in- 
stantaneous potential can be consistently derived from the underlying field theory 
of QCD. 

2. Defining the static heavy quark potential 

What makes heavy quarkonia particularly suitable to theoretical treatment is the 
inherent separation of scales between its constituent mas^^ {rric — 1.29GeV, 777,5 — 
4.6GeV) and other typical scales in relativistic heavy ion-collisions. Compared to 
e.g. the temperature of the deconfinement transition Tc ^ 200MeV, the intrinsic 
scale of QCD Aqcd ^ 200MeV and the typical momentum exchange in the bound 
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state p, we find 

(1) 

rriQ rriQ tuq 

This in turn tells us that neither thermal nor quantum fluctuations can sponta- 
neously pair create a heavy Q and Q and thus a non-relativistic description can be 
attempted. 

Whenever there exists a hierarchy of scales, effective field theory (EFT) provides 
a general framework on how to reduce the complexity of the description by treating 
explicitly only the physics at the energy scales of interest. (For a review of the 
technique applied to QCD at T=0 see Refs.'^-'^'^^'and references therein.) In our case 
the physics of the bound state is well separated from the hard scale mg, which 
we hence wish to integrate out. The notion of "integrating out" encompasses three 
actions: First we need to identify the relevant degrees of freedom at the energy 
scale of interest; secondly the most general Lagrangian in these new variables is to 
be constructed while retaining the symmetries of the underlying field theory. Last 
but not least, the residual influence of the higher lying energy scale is incorporated 
through the process of matching, where suitable correlation functions in the original 
and effective language are made to agree at a certain energy scale. 

2.1. An effective field theory of Pauli spinors 

The first step on the path towards a heavy quark potential thus lies in applying the 
above concepts to QCD. We begin by integrating out the hard scale, which leads 
us to the effective field theory NRQCD. The starting point is the QCD Lagrangian 

>^QCD = 

- ^FrK^ + E + WA,) - m[)q' ^Q{ir{d, + igA,) - mg)Q, 

written in terms of matrix valued gauge fields A^ = A^^T^ ^ light medium quarks q 
and heavy quark fields denoted by Q. 

As the physics below the hard scale does not involve pair creation of heavy 
quarks, the upper and lower components of the Dirac spinor Q = (^, x) themselves 
suffice to capture the relevant physics and we proceed to separate ^ and x by means 
of a Foldy-Tani-Wouthuysen (FTW) transformatioriMH The resulting NRQCD 
Lagrangian presents itself as a series of terms of increasing powers in the inverse 
heavy quark mass tUq^^ the first few of which read^^ 

+ X^(zI)o-mg + ^D2 + ...)tx. (2) 

Note that in addition to the transformation of fields, we have inserted so called 
Wilson coefficients q in each term that scales explicitly with a power of m^^. These 
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complex numbers encode the remnants of the physics at higher energy scales, such as 
e.g. a possible contribution of a gluonic cloud to a kinetic mass of the quarks via ci. 
One difference to the usual NRQCD approach is that here we do not integrate out 
gluons and light quarks with hard momenta, as we assumed Tmed <^ and our aim 
is to explicitly incorporate their contribution through the use of non-perturbative 
lattice QCD simulations. For our purpose of deriving a static and spin-independent 
potential, we can safely set all values of q to unity in the following. (For further 
information on the explicit determination of the q's through matching we refer the 
reader to Ref.l^. 

The strength of the FTW transformation is to yield a systematic expansion 
in the inverse rest mass, which at the same time alerts us to the limits of any 
potential picture derived via NRQCD. Since after summing up an infinite number 
of terms, the underlying QCD Lagrangian emerges, couplings between the upper 
and lower components x ^ will necessarily reappear at some order and the 
absence of explicit pair creation is not guaranteed. The first relevant terms, such as 
i^XX^i indeed appear at order ^rig^, signaling the break down of the non-relativist ic 
approximation at higher orders. 

The notion of potential has so far not found entry into the description of the 
dynamics, as the Pauli spinor fields ^ and x interact via the explicit mediation of 
the gauge fields A^. One possibility to introduce V{r) is to integrate out all degrees 
of freedom down to the so called soft scale, characterized by energies E ~ rnqv. The 
Wilson coefficients in this approach are non-local and represent potentials between 
well defined color singlet and octet configurations of the heavy quarks. Termed 
pNRQCD, this effective field theory allows a consistent power counting and the 
choice of degrees of freedom in the form of correlated pairs of heavy quarks has 
proven successful at zero temperature^^ as well as in in-medium applications^^. 

On the other hand, in the quark-gluon plasma we expect the heavy quark bound 
states to melt eventually and thus ultimately our description should comprise as 
relevant degrees of freedom also decorrelated pairs of individual heavy quarks. 



2.2. Quantum mechanical path integrals 

The question remains how to describe the heavy quark bound state while empha- 
sizing the propagation of its individual constituents? One possible way is to turn 
to non-relativistic path integral^IEU, which operate with point particles^^ along 
fiuctuating trajectories and their conjugate momenta p^. In particular they will 
allow us to to read off the interaction potential by applying the transfer matrix 
prescription to the complex weighting factor exp[z J dt{pz — H)] of a two-body path 
integral. (For an application to a single heavy quark see Ref.^^) 

Starting point is a suitable correlation function in field theory that serves as 
propagator for the heavy quarkonium system. Here we choose the forward correlator 

^ as it encodes how probable it is to find a bare QQ state at time t after starting 
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the evolution at time t=0 
^^(xi,yi,x2,y2,t) 



(3) 



(Q(x2,t)[/(x2,y2,t)Q(y2,t)Q(yi,0)/7t(xi,yi,0)Q(xi,0)) 



For the expression to be gauge invariant we connected the point split fermionic 
fields by straight spatial Wilson lines [/(x, y,t). 

The results of Sec 2.1 allow us to to reexpress Eq.(|3| in the Pauli spinor fields 
of the NRQCD effective field theory 

^NRQCD(xi,yi,X2,y2,t) = 

^K, ^]D[x. X]D[A, q, g]^t(x2, t)/7x(y2, t)x^(yi, 0)/7t^(xi, 0)e^^--Q^- . (4) 



The trick here is to actually carry out analytically the quadratic integration over the 
heavy GraBmann fields, which leaves us with non relativistic propagators, describing 
the evolution of quark and anti-quark in the background of the medium gluon field 
A^. As heavy quarks in NRQCD, the fields ^ and x cannot participate in virtual 
loops and thus no fermion determinant has to be taken care of. In practice we 
have to replace pairs of and XX^ by the Green's functions K and ^ i.e. the 
functional inverses of the NRQCD action defined by 



idt - gAo 



lim K(xi,X2)=(5('Hxi-X2). 

Xi^X2 

It is at this point that the connection to quantum mechanical path integrals is made, 
as the propagator of a particle with trajectory z and momentum p can be expressed 



i^(xi,X2) = J Vz J Pp T exp i J ds (^pz 



2mQ 



exp 



ig / dy^A^iy) 



Since quark and anti-quark each enter through a separate path integration, the 
forward propagator in the language of quantum mechanics takes the form 



-2imQt 



ry^2 ry2 
/ P[zi,pi] / P[Z2,P2] 
Jxi Jyi 



(5) 



Texp 



dx ^A^i^x) 



where we have separated a field independent term containing the kinetic terms, 
quadratic in the momenta, and a field dependent term, which still requires us to 
average over the medium degrees of freedom. 

In order to understand how a potential of the two-body system emerges from 
Eq.([5|, we focus on its last term, which is nothing but the thermal real-time Wil- 
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son loojQ If a potential picture were applicable at all times, we expect to be 
able to rewrite the expression as an exponential over a time independent func- 
tion V{r = |z2 — zi|). Since what we do however amounts to replacing a retarded 
field theoretical interaction between quark and anti-quark by an instantaneous non- 
relativistic potential, we have to take into account the possibility that the function 
F(r, s) does actually vary at early times before approaching a constant value only 
at late times 



Wc\{r,t) = (rexp ig J dx^A^{x) ^ = exp -i J dsV{r,s) 



(6) 



Note that we assume that the heavy mass limit allows us to neglect terms connect- 
ing different times along the path of the quark-antiquark pair in Eq.([6|. Going to 
late times finally leads us to the following defining equation, connecting the time 
evolution of the real-time Wilson loop to the static inter-quark potential 

V{r) hm . (7) 

^^oo Wa{r,t) 

3. Extracting the static heavy quark potential 

Now that the field theory of QCD in form of the real-time Wilson loop has been 
connected to a non-relativistic potential description of heavy quarkonia in Eq.Q, 
we have to ask how the actual values of such a potential can be determined in 
practice. 

One possibility, worked out in detail in the ground-breaking contribution by 
Laine et. al.^^ is to use resummed perturbation theory, the so called hard thermal 
loop approximation^^*H^ in describing the medium the heavy quarks travel in. This 
gauge invariant prescription allows us to sum an infinite number of Feynman dia- 
grams for thermal gluons and light quarks and already captures essential features 
of the in- medium physics at high temperatures. Evaluating the Wilson loop at first 
non-trivial order in the gauge coupling and inserting into Eq.Q yields^^ 

^HTL r = rriD + - m^r , 8 



Stt 

= 2 / dz 



^ ^ " sinlzx] 







1 



zx 



{z^ + 1)2 

One finds that the potential is complex with the real part exhibiting a Debye 
screened form with screening mass mjj = ^^T^^^ + The appearance of 

an imaginary part on the other hand can be linked to collisions between the light 
partons in the medium and the gluon mediating the interaction between the Q and 
Q, i.e. the phenomenon of Landau damping^^. While Im[yHTL] grows quadratically 

^Even though this quantity is usuahy referred to as thermal, it is not periodic in imaginary time. 
Intuitively this is a result of the static quarks not being thermalized. The gluons of the thermal 
medium and their periodicity however introduce an upward trend into the Euclidean time data, 
seen most clearly as r — /3. 
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at small m^^r < 1 it saturates to a constant value at large m^^r 1, which in 
turn can be interpreted as twice the energy loss of a single quark traversing a heat 
bath. First corrections to Eq.([8| have been presented in a pNRQCD contextP^, in- 
corporating the first oder of the multipole expansion. There it was shown that the 
breakup of a singlet to an octet state under the influence of a color electric field can 
also contribute to the imaginary part. 

While asymptotic freedom guarantees the correctness of weak-coupling results 
at high temperatures, we expect the quark-gluon plasma to be strongly interacting 
in the phenomenologically important region around the phase transition^^. One 
novel possibility to, at least qualitatively, elucidate such non-perturbative QGP 
physics is the gauge-gravity dualit}'^^. By mapping a certain type of strongly coupled 
conformal Yang-Mills theory to a weakly coupled classical field theory in higher 
dimensions, calculations that are technically impossible in the former can actually 
be carried out in the latter. This approach has recently been used^^ to determine the 
values of the thermal real-time Wilson loop and subsequently the static potential 
according to the definition Eq.Q. While the obtained real part at small distances 
shows Coulombic behavior, at large distances it exhibits a rather long tail, whose 
value decreases with increasing temperature. The imaginary part on the other hand 
while being zero up to a threshold rth = 0.62 (ttT)"^ runs linearly in distance already 
dominating the real part at relatively short rdom = 1.72(7rT)~^. 

In the following we investigate the non-perturbative aspects of heavy- 
quarkonium physics, while staying firmly within the framework of QCD. One pos- 
sible way to do so is to rely on Monte Carlo simulations of spatially regularized, 
so called lattice QCD. This approach allows us to evaluate observables at any tem- 
perature, with the only but significant limitation that all calculations have to be 
carried out in imaginary time instead of real-time. Hence no direct access to the 
real-time Wilson loop of Eq.Q is possible. In order to circumvent this obstacle and 
to obtain the sought after information necessary to evaluate the potential, we will 
introduce a spectral representation of the Wilson loop insteadHEH 



3.1. A spectral representation 

It has been showrP^ that the Wilson loop admits a spectral representation, which 
amounts to nothing but a Fourier transform over a positive definite function pn(r, uj) 

/CO 
due-"^'pair,w). (9) 
-OO 

The benefit of writing the Wilson loop in this way is that its time dependence only 
appears in the kernel of the integral. If the spectral information is known in the 
form of pn(r, cj), changing from real- to imaginary time just means to go over from 
a Fourier to a Laplace transformation 

pOO 

Wa{r,T)= / doje-'^^pa{r,uj). (10) 
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Technically there is however an important difference, since an exponentially damp- 
ing Laplace kernel makes the inversion of the relation between spectrum and Wilson 
loop much more demanding than an oscillating Fourier kernel. And it is exactly such 
an inversion that we have to perform if we wish to extract from the noisy lattice 
QCD estimates of VKn(r, r) the spectral information to reconstruct VKn(r, (See 



Sec. 3.2). 



The question to ask here is whether we can already extract the values of the 
potential from spectral information of the Wilson loop? To find out we follo^^^ by 
inserting the spectral representation of Eq.([9| into the defining Eq.([7|) and obtain 

V{r) = lim \ (11) 

^^oo duj e-*^^ (r, uj) 

Even though we know from Eq.Q that the potential is only connected to late time 
physics and naively speaking Eq.(pT]) hence signals that we will need to consider 
spectral features only at small frequencies, the fact that the Fourier transform relies 
on late and early time information needs to be kept in mind. 



Indeed from our derivation of the potential in Sec 2.2 we know that the de- 



scription of the time evolution of the Wilson loop in the form 

idtWn{r,t) = <i>{r,t)Wn{r,t) (12) 

in general requires a complex function <l>(r, t). Only at late times it approaches a 
constant, which we identify with the potential 

lim <l>(r,t) = V{r). (13) 

Following Ref.'^ we introduce an additional function (j)(r^t)^ which encodes the 



deviation of <l>(r, t) = V{r) + (l){r^t) from the potential and solve Eq.(12) to obtain 
a general expression for the Wilson loop 



Wu{r,t) exp 



■ i {Re[V] {r)t + Re[cr] (r, t)^ - \lm[V] {r)\t + Im[cr] (r, t) 



The quantity <j(r, t) = (j){r^ t)dt and its asymptotic values are defined as croo{r) = 
a{r^ \t\ > tgo) = (j){r^t)dt. Carrying out the Fourier transform^^ tells us that 
the most general spectral shape we will encounter at low frequencies in pn(r, cj) 
reads 

n (. , .^ - 1 imk^Kr) |Im[y] (r) |cos[Re[aoo] (r)] - {Re[V\{r) - uj)sm[Re[cj 



TT Im[y](r)2 + (Re[y](r) - Lof 

+«o(r) + «i(r)tQQ(Re[y](r) - lo) + «2(r)i|Q(Re[y](r) - ujf + ■ ■ ■ (14) 

The first term in the above Eq.([T4| takes the form of a skewed Lorentzian embedded 
in a polynomial background, characterized by real numbers ni{r). Even though 



^Note that even though the Euchdean time quantities are all purely real, the presence of the 
Fourier kernel makes it possible for the potential to have an imaginary part in the end. 
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the result in Eq.(14) appears quite involved, inserting it into our defining formula 
Eq.(ll) and carrying out the contour integration in the lower half of the complex 
plane yields a time independent potential with V{r) = Re[V](r) — z|Im[V](r)|. In 
earlier studied the spectral features were assumed to be a perfect Breit-Wigner, 
which we can now understand from the above derivation as corresponding to the 
additional assumption dt^{r^t) = , i.e. that the potential picture be applicable at 
all times. 



3.2. The maximum entropy method 

All conceptual ingredients are now in place to connect the spectral information of 
the Wilson loop to the concept of a static inter-quark potential. One remaining issue 
of practical importance is how the spectrum is actually measured in lattice QCD. 
Over the last decade a Bayesian approach to this question has found widespread 
adoption, the Maximum Entropy MethocHZEH (MEM). What we are facing in the 
context of lattice QCD when attempting to invert Eq.([To| is a inherently ill-defined 
problem. The reason is that the Monte-Carlo estimates of VFn(r, r) are available only 
on a discrete number of points and carry a finite uncertainty. If we were to attempt 
an extraction of a continuous function pn(r, cj) via fitting, we immediately find 
that an infinite number of degenerate solutions will ensue, all of which fit the data 
within their errorbars. 

Based on Bayes theorem, we can nevertheless give meaning to the inversion 
task by emphasizing the role of prior information /. The probability of a test spec- 
tral function to be the correct spectral function, given measured data D and prior 
information /, can be written as the product of two terms 

P[p\D,I]o:P[D\p]P[p\I]. (15) 

The first one denotes the so called Likelihood probability, the second one refers to 
the prior probability. While P[D\p] is nothing but the usual fitting term, the 
MEM uses the so called Shannon- Jaynes entropy S in P[p\I] to incorporate prior 
information, in the form of a function m{uj)^ while enforcing the positive definiteness 
of p 

Pmem[p\I{^)] ^ exp[a5] = expU^ (pi - mi- p/log[^] j . (16) 

L \ mi 

If we set out to determine the most probable spectral function, we are lead to 
numerically solve the following stationarity condition 



^P[p\D,I{m) 
opi 



cx ^(^P[D\p]Pmem[p,I]) 

P=Pmem 



= 0. (17) 

P=Pmem 



Here the otherwise underdetermined fitting, which attempts to solely reproduce 
the measured data, competes with the prior term which tries to bring the spectrum 
as close to the prior function m{uj) as possible. It can be shown^ that Eq.(17) 



actually possesses a unique solution Pmem if such a solution exists. Intuitively the 
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Fig. 1. A schematic overview of the strategy to non-perturbatively extract the static heavy- 
quark potential from lattice QCD simulations. From measurements of the Euclidean Wilson loop 
at different separation distances along the imaginary time axis, we can obtain the corresponding 
spectral information by applying the Maximum Entropy Method. Based on the most general 
functional form given in Eq.( |14| > the lowest lying spectral feature at positive frequencies is fitted 
and the values of Re[y](r) and Im[y](r) are read off. 



reason is that by taking together measured datapoints and the supphed prior func- 
tion, we have at our disposal more points of information than parameters to extract. 
Note that in the extreme case of no measured datapoints, by definition, the function 
771(0;) itself represents the correct spectrum. 

Consequently, parts of the spectrum Pmem are constrained by the measured 
data, parts of it by the provided prior function. The introduction of such a regular- 
izing mechanism into the fitting hence allows us to unambiguously identify what 
spectral features are reliably encoded in the measured data, i.e. by repeating the 
MEM extraction with many different functional forms of the prior m{uo). For tech- 
nical details on the numerical implementation and the question of an appropriate 



search space in which the extremum, defined by Eq.(17), can be located, we refer 
the reader to RefsEMIMEIl 

It is here that we can finally put all pieces of the puzzle together, as sketched in 
Figjl] and present the current state of knowledge on the real and imaginary part of 
the heavy-quark potential around the phase transition from quenched lattice QCD 
(Fig {2]). The shown values are based on measurements of the Euclidean Wilson 
loop^ in a purely gluonic medium on anisotropic lattices of size 20^ x 12, 24, 36 at 
a bare anisotropy ^5 = 3.2108. The choice of /3 = 6.1 {ax = 0.097fm) corresponds 
to the temperatures T = 2.33Tc, l.lTTc and 0.78Tc. After carrying out the MEM 
on the imaginary time data at different separation distances, the lowest lying pos- 
itive peak is fitted by the shape introduced in Eq.(14) (see also Ref.^^) and the 



corresponding values for Re[]/](r) and Im[]/](r) are read off. 

In the hadronic phase, the real part of the potential Re[Fn] agrees with the 
color singlet free energies F^{r) within its errorbars. Note that this is not obvious, 
since the free energies are obtained from the correlator of Wilson line^ W|| (r, r) in 
Coulomb gauge at a single point in time r = /3. In the spectral picture (see Figjl]) 

^This quantity ensues if the straight spatial Wilson lines ?7(x,y,t) are removed from Eq.(|3j 
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Fig. 2. The real and imaginary part of the static inter-quark potential in quenched lattice QCD, 
based on the MEM spectra of Euclidean Wilson loop measurements. We fit the spectrum with 
Eq.(|14|, taking into account skewing and up to quadratic background terms (SQL). While as 
previously indicated, the real part at T = 0.78Tc coincides with the color singlet free energies 
F-^(r), the error bars are still too large to reach a definite conclusion on whether or how the real 
part at T > Tc exhibits Debye screening. Note that the artificially strong rise from the Breit- 
Wigner fit disappeared at T = 2.33Tc. Finite values of ImfT/g] below Tc are most probably due 
to a finite resolution of the MEM, introduced by statistical uncertainty in the data. In the QGP 
phase however the curvature of W\j{r,T) leads to finite values. 
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Fig. 3. The real and imaginary part of the static inter-quark potential in quenched lattice QCD, 
based on the MEM spectra of Euclidean Wilson line measurements in Coulomb gauge. We fit the 
spectrum with Eq.(|14|, taking into account skewing and up to quadratic background terms (SQL). 
A good signal to noise ratio allows us to see both the linearly rising nature of the real part below 
Tc as well as a clear sign of screening in the values above 2Tc- The apparent finiteness of the 
imaginary part below Tc needs to be attributed to a finite resolution of the MEM, introduced by 
statistical uncertainty in the data. In the QGP phase at T = 2.33Tc the finite imaginary part 
seems to show a relatively strong rise at r < 0.4fm while fiattening off at larger distances. 



this point is not connected to the potential peak at positive frequencies but instead 
to an exponent iahy suppressed feature at negative frequencies. 

Above the deconfinement transition the large error bars on the reconstructed 
values prohibit us from making a conclusive statement on whether or not the po- 
tential shows Debye screening, as predicted from HTL calculations. We find that 
with the the fitting function Eq.(14) introduced in Ref.*^, the counterintuitively 
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large rise in the real part observed in Ref.^^ at T = 2.33Tc is significantly reduced. 
What is important is that for the imaginary part Im[]/n] we find finite values in the 
deconfined phase, even after subtracting the values found below Tc, as a baseline for 
the limited resolution of the MEM. Inspecting the Euclidean time data of VFn(r, r) 
at T > Tc, as a cross-check, reveals a significant curvature at intermediate values 
of r which then translates into a width in the spectrum. This gives us confidence 
that with these extracted values at T = 2.33Tc the existence of an the imaginary 
part in the quark gluon plasma is confirmed on the lattice. 

If we interpret the Wilson loop as the correlator of two temporal Wilson lines 
in a spatial axial gauge, we might attempt to redo the extraction of the potential 
also in Coulomb gauge based on W||(r, r). This observable exhibits a much higher 
signal to noise ratio on the lattice. As shown in Fig|3] it thus allows for a much 
more accurate determination of the spectral information, which in turn leads to 
smaller errorbars in both Re[l^|] and Im[V||]. The real part appears to move from a 
confining linear rise, similar to F^(r), below Tc, to a screened form, which however 
does not coincide with the free energies anymore. The finite values of the imaginary 
part in the hadronic phase are still to be attributed to an artificial resolution limit 
in the MEM, coming from the statistical uncertainties of the Euclidean data. The 
imaginary part at T = 2.33Tc on the other hand, which is related to a signal 
encoded in the Euclidean correlator, data shows a rise at small distances r < 0.4fm 
before an apparent fiattening sets in at larger distances. 

The presented results for the static inter-quark potential from lattice QCD tell us 
that while conceptually a clear path has been established connecting the quantum- 
chromodynamics of a heavy quark and antiquark in a thermal medium with a po- 
tential picture, we need to improve the reliability of the extracted values. One route 
which needs to be taken is to measure the Wilson loop to much higher accuracy 
on even larger lattices. As the current calculations are still based on the quenched 
approximation the multilevel algorithnP^ might be the appropriate technical tool to 
use. The MEM itself is known to introduce an artificial width into the reconstructed 
spectra as its resolution is limited by the amount of uncertainty in the supplied data 
sets. Hence it needs to be established through mock data analysis how precise we 
need to measure l^n(r, r) in order to reliably reconstruct the spectral width even 
below the deconfinement transition. Last but not least, we should extend the strat- 
egy reviewed here to measurements of the Wilson loop in the context of dynamical 
QCD, where significant effects of deconfined light quarks to the screening of the real 
part and strength of the imaginary part are expected. 

4. Heavy quarkonium as open quantum system 

The previous sections have shown that the static inter-quark potential in the de- 
confined phase most certainly exhibits an imaginary part. We might ask how such 
an apparent non-hermiticity can be interpreted in the context of the in-medium 
evolution of the heavy-quark bound state. To this end let us remember that we 
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based the definition of ttie potential on tiie forward correlator in Eq.(|3| and 
the presence of a finite Im[V] tells us that it is these correlations, which damp away 
with time. This however does not mean that the constituents of the heavy quarko- 
nium themselves disappear, which they cannot, since annihilation into hard gluons 
is not taken into account when operating to order ttIq^ in Eq.([2|. While the overall 
strength of the correlations is reduced by Debye screening in the presence of decon- 
fined partons, it is scattering with the light quarks and gluons, which actually leads 
to a decrease in correlations over time. Such a loss of correlation in turn manifests 
itself as an imaginary part in the Schroedinger equation for . This goes so far 
that after some time, changes in one particle do not infiuence the state of the other. 
At this point the notion of a bound state becomes devoid of meaning, the heavy 
quarkonium has melted. 

To formalize this idea of decoherence in the language of quantum mechanics and 
to see how the imaginary part arises from the thermal fiuctuations in the medium 
surrounding the QQ, we turn to a description based on the theory of open quantum 
systems. (Recent work in this direction can be found in Kefs. |43|44|45|46|47|) ^j^-g ^^^^ 
established framework provides the conceptual tools to describe the infiuence of 
a medium onto a small subsystem, a topic thoroughly investigated in condensed 
matter theor}/^^. Assume that both the medium and the QQ can be described 
quantum mechanically, so that the overall Hamiltonian i^fuii = ^/uii is hermitian 

d 
It 

i.e. states evolve unitarily and the density matrix of states (j{t) follows the von 
Neumann equation. If we now wish to describe the system solely in terms of the 
heavy Q and Q, as we have done in our attempt to derive an effective Schrodinger 
equation, we have to trace out all other degrees of freedom in the system. Their 
infiuence on the evolution of the subsystem manifests itself in the appearance of 
a stochastic element, such as noise, both in the master equation of the density 
matri:jP^ as well as in the evolution equation of the wavefunction. 

Since a Schrodinger equation does not possess a notion of thermal fiuctuations, 
one necessarily goes over to an ensemble of wavefunctions ^ qq so that the corre- 
sponding density matrix of states in the subsystem can be expressed as 

aQQ(t,r,r') =Twd[^(i,r,r')] = (^^^(r, t)vl/^Q(r', t)). (19) 

By definition, decoherence in this context represents the phenomenon that the in- 
teractions with the surroundings select a certain basis of states in the QQ system 
in which the density matrix cfqq becomes diagonal after the decoherence time tdc 
has passed. To make as close contact as possible with the potential extracted from 
lattice QCD, we will however study the infiuence of the interaction with the thermal 
medium directly on the level of the wavefunction. 

In Sec |2.2| we connected the spectral features of the Wilson loop to a complex 
potential. Here we take a different route^as indicated in FigJZJ The intuitive idea is 



HiuW = Hmed ^ Iqq + ^med ^ HqQ + i^int, :^Cr(^) = -^[^full, Cr(^)], (18) 
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Effective Field Theory Open Quantum System 



Fig. 4. (left) Reinterpreting the spectral features in the Wilson loop spectrum. Instead of as- 
signing the width to an imaginary part of the potential as done in the EFT approach, the open 
quantum systems approach works with a purely real potential which is perturbed by the thermal 
medium. It is the strength of the thermal noise r(r, r) which is then characterized by the width, 
(right) Schematic view of the difference between a description on the level of the full system with 
a hermitian Hamiltonian H = Hqq + ii^med + -f^int cind a stochastic potential description of the 
QQ system after tracing out the medium degrees of freedom. 



that the thermal fluctuations, i.e. light quarks and gluons in the QGP, will perturb 
the potential acting between the heavy Q and Q at each step in the time evolution. 
The average of this purely real potential VQQir) corresponds to what was previously 
called Re[y], while its variance is to be connected to the width of the spectral 
features. Based on this paradigm, let us set out to construct a fully unitary time 
evolution operatoil^ for each microscopic realization of the wavefunction ^ qq in 
the ensemble 



^QQ(r,^) = Texp 



/ d5 I - — + 2mg + Vqq{v) + e(r, s)] 
Jo ^ ^ 



*Qc3(r,0). 



(20) 

For the sake of simplicity we restrict ourselves to a Markovian noise term 
(0(r, t)) = 0, which however carries a non-trivial spatial correlation structure 
(0(r, t)0(r', t')) = ^(^t^t/r(r, r') characterizing the thermal properties of the QGP 
medium. An equation of motion for the heavy quarkonium is obtained once the 
operator of Eq.(20) is expanded following the rules of stochastic differential Ito 
calculus 

z^^gg(r, t) = ( - ^ + 2mg + Fgg(r) + e(r, t) - i^e\r, t)) ^gg(r, t). (21) 

At this point we are in possession of a stochastic Schrodinger equation that pre- 
serves the norm of each individual wavefunction of the ensemble while introducing 
decoherence in the physical heavy quarkonium state, as seen from taking the average 

4(*QQ(r,t)>= (-^ + 2mQ + FQQ(r)-^r(r,r))(*QQ(r,t)). (22) 

Thus within this open quantum systems approach, the imaginary part of the EFT 
approach emerges naturally from the diagonal correlations of the thermal fluctua- 
tions. On the other hand, extracting the imaginary part of the EFT potential from 
lattice QCD allows us to pinpoint parts of the noise structure of the QGP. 
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To describe the suppression of heavy quarkonia, we have to state first what 
it is that we measure in experiment. Naively speaking, when heavy quarkonium 
is created in the partonic stage of a cohision, it appears as eigenstate ^n(r) of a 
vacuum Hamilton operator H^q. After entering the QGP, the stochastic evolution 
will lead to a reshuffling of states but the particle we finally observe in the detector 
is still a vacuum eigenstate. Hence we wish to ask how probable it is to find such a 
vacuum state after a certain time of evolution in the thermal medium, a question, 
which is answered by the following projection 



(t) = / dVdV rnir){^QQ{r,mQQ{r',t))4>n{r'). (23) 



Since this quantity is connected to the density matrix of states introduced in Eq.(|19[) 
it will be susceptible also to the off-diagonal elements in r(r,r') in contrast to the 



averaged wavefunction of Eq.(22). By observing the time evolution of Cnn{t) we are 
confident to learn, within the limitations of the non-relativistic approach, whether 
and how a heavy quark bound state melts as it evolves in a thermal medium. 
Conversely by comparison with experimental suppression data, one might attempt 
to infer properties of the QGP itself. 

Obviously we are only at the beginning of a thorough understanding of the 
open quantum systems nature of heavy quarkonium. One immediate question is 
e.g. how to extract not only the diagonal but also the off-diagonal components of 
r(r, r') from lattice QCD. A more fundamental challenge to the presented approach 
is the fact that we use a static potential to describe heavy but finite mass quarks in 
Eq.(20). As was recently shown in a perturbative rederivation and generalization^ 



of the real-time dynamics inspired by the Feynman Vernon infiuence functionaP^, 
we are missing a crucial element. In fact, dissipation, without which the heavy 
quarks can never thermalize, must be included to obtain e.g. consistent Ehrenfest 
relations. Once thermalization has been achieved it will be of interest to compare 
the resulting distribution of states to the spectral functions from HTL perturbation 

theor^EEMlland lattice QCa54|55|56|39|57|58|59|6m 

In addition, a possible route to include explicit color degrees of freedom has been 
presented^, a feature we require to allow the heavy quarkonium to melt consistently 
into its constituents. Further studies in this direction should be undertaken as they 
open the way towards a consistent description of both single heavy quarks and their 
bound states in the QGP. 
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